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Abstract We present self-similar solutions for radiation inefficiently accretion flows 
(RIAF) around black holes in the presence of outflow and global magnetic field. The 
influence of outflow is taken into account by adopting a radius dependent of mass accre- 
tion rate M = M^(r/r^) s with s > 0. Also we consider convection through a mixing 
length formalism to calculate convection parameter a con . Our numerical results show that 
by increasing all components of magnetic field , the surface density and rotational veloc- 
ity increase although the sound speed and radial infall velocity of the disc decrease. Also 
we have found out that the existence of wind will lead to reduction of surface density as 
well as rotational velocity. Moreover the radial velocity, sound speed, advection parame- 
ter and the vertical thickness of the disc will increase when outflow becomes important in 
the radiation inefficiently accretion flow. 
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1 INTRODUCTION 

Observations of accreting black holes of different mass show impressive similarities of data which point 
at identical physical process in the accretion flow. Also, black hole accretion discs demonstrate a great 
variety of physical conditions, so we may have a variety of accretion regimes. Existing theories describe 
different regimes of black hole accretion flows, which can be realized under different physical condi- 
tions. Accreting black holes in nearby galactic nuclei and low-state X-ray binaries are much dimmer 
than the standard Shakura-Sunyaev disc model would predict. A phenomenon of under-luminous ac- 
creting black holes in X-ray binaries and super-massive black holes in galactic nuclei has stimulated 
the recent investigations of radiation inefficient accretion flows (RIAFs) (see Narayan, Mahadevan & 
Quataert 1998, Kato, Mineshige & Fukue 2008 for reviews). In such a flow, radiative losses are small 
because of low particle density of accreting flow at low accretion rates. Contrary to standard Shakura- 
Sanyaev disc model, which successfully explains radiation soft and luminous X-ray sources, models of 
RIAFs are used to explain significant deficit of radiation observed in some X-ray sources. A particular 
example of such under-luminous sources is the Galactic center, Sagittarius A*, with host a 2 x 10 6 solar 
mass black hole. The Galactic center has a luminosity that is well bellow the estimated value based on 
Shakura-Sanyaev accretion disc Model (Melia & Falcke 2001). 

ADAFs have an opposite regime to that of the standard model. In the standard model, the flow is 
described in such a way that the heat generated by the viscosity radiates out of the system immediately 
after its generation (Shakura & Sunyaev 1973). These advection-dominated accretion flows occur in two 
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regimes depending on their mass accretion and optical depth. Actually, the optical depths of accretion 
flows are highly dependent on their accretion rates. In a high mass-accretion rate, the optical depth 
becomes very high and the radiation generated by the accretion flow can be trapped within the disc. 
In this case, the optical depth is very large, and photons, which carry the most of internal energy, are 
trapped inside the inflowing matter and can not be radiated away. This type of model is called 'slim 
disc', or optically thick ADAF. Although the radiative efficiency of optically thick ADAF is also low, 
we usually only call the optically thin ADAF as RIAF. In the limit of low mass-accretion rate, the disc 
becomes optically thin. In this case, the cooling time of accretion flows is longer than the accreting 
time-scale. The energy generated by accretion flows therefore remains mostly in the discs, and the 
discs cannot radiate their energy efficiently. This type of accretion flow is named a radiation-inefficient 
accretion flow (RIAF). This type of accretion flow has been investigated by many authors (Narayan & 
Yi 1994; Abramowicz et al. 1995; Chen 1995). At the same time as the ADAF model was proposed, 
it was realized the ADAFs are likely to be convectively unstable in the radial direction because of the 
inward increase of the entropy of accreting gas (Igumenshchev & Abramowicz (1999), Stone, Pringle & 
Begelman (1999), Igumenshchev, Narayan & Abramowicz (2003)). Most recent work focusing on the 
convective instability is by Yuan & Bu (2010). Two and three dimensional simulations of low-viscosity 
RIAF have confirmed the convective instability in these flow (Igumenshchev, Chen & Abramowicz 
1996; Igumenshchev & Abramowicz 2000, McKinney & Gammie 2000). Narayan, Igumenshechev & 
Abramowicz (2000) and Quataert & Gruzinov (2000) construct another analytical model of RIAFs 
which was based on a self-similar solution which was called the convection dominated accretion flow 
(CDAF). CDAF consist of a hot plasma at about virial temperature and have a flattened time-averaged 
radial density profile, p oc r -1 / 2 , where p is the density and R is the radius. In CDAFs the most of energy 
which released in the innermost region of accretion flow is transport outward by convection motion. 

The self-similar CDAF model as the same as other self-similar models is very clear and instructive, 
but it has some limitations. It is only a local, not a global solution for RIAF, it the sense that it can only 
be valid for a region far from boundaries. So it can not reproduce physical behavior of accretion flow 
it it's transonic radial motion- the most fundamental feature of black hole accretion flow. Abramowicz 
et al (2002) did suggest two-zone structure for RIAF: an outer convection dominated zone and an inner 
advection dominated zone separate at a transition radius ~ 50r g . 

There is some observational evidence that accretion process is often and perhaps always associ- 
ated with wind. Mass loss or wind appears to be a common phenomenon in many astrophysical sys- 
tems. These wind and outflow mechanisms are observed in micro-quasars, YSOs ( Ferrari 1998, Bally, 
Reipurth & Davis 2007; Whelan et al 2005). It is belief that disc wind/outflow contributes to loss of mass, 
angular momentum, and thermal energy from accretion discs (e.g. Piran 1977; Blandford & Payne 1982; 
Foschini 201 1, Knigge 1999). Various driving sources are proposed, such as thermal, radiative and mag- 
netic mechanism. The name of the wind depends on its driving mechanism. In this paper we will follow 
the hydrodynamical (thermal) wind which has been discussed by many authors (e.g. Meier 1979, 1982; 
Fukue 1989; Takahara, Rosner & Kusnose 1989). 

In this paper we will discuss the properties of CDAFs in a general large-scale magnetic field with 
hydrodynamical wind. We will concentrate on the self-similar solution. This study is motivated by re- 
cent works of Zhang & Dai 2008 who showed that the effect of large scale magnetic field on the CDAFs 
with a constructive self-similar solution. Because of CDAFs usually modeled for outer regions of RIAFs 
where the hydrodynamical outflow becomes important, we investigate the role of outflow in the dynam- 
ical structure of magnetized CDAFs. We will show the basic equations and self-similar solution in next 
sections. 

2 BASIC EQUATIONS 

We are interested in analyzing the structure of a magnetized hot accretion flows bathed in a global 
magnetic field where convection and wind play an important role in energy and angular momentum 
transportation. So we suppose a rotating and accreting disc around a compact Schwarzchild black hole 
of mass M*. Thus, for a steady axi-symmetric accretion flow, i.e., d/dt = d/dcp = 0, we can write the 
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standard equations in the cylindrical coordinates (r, ip, z) centered on the accreting object. We vertically 
integrated the flow equations and, all the physical variables become only function of the radial distance 
r. Moreover, we consider a magnetic field in the disc with three components, B r , B v and B z . We also 
neglect the relativistic effects and Newtonian gravity in radial direction is considered. 
The equation of continuity will be 

— r£u r ) + - — =0 (1) 

Or 2n Or 

where £ is the surface density at the cylindrical radius r, which is define as £ = 2pH, p being midplane 
density and H the disc half-thickness and v r is the radial infall velocity. Also the mass loss rate by 
outflow/wind is represented by M w , so 

M w (r) = I 4irr'm w (r')dr' (2) 



where m w (r) is the mass loss rate per unit area from each disc face. We can write the dependence of 
accretion rate as follows (Blandford & Begelman 1999; Shadmehri 2008) 

M = -2-KrEvr = Mo( — Y (3) 

where Mq is the mass accretion rate at the outer radius of the disc ro and s with s being constant of 
order unity (Blandford & Begelman 1999). 

Considering equations ([TJ-fO, we can write 



M , r , , 

m w = s-^(-y- 2 (4) 



The equation of motion in the radial direction is 

where v v is the rotational velocity, c s is the sound speed, which is define as (? s = p gas / p, Pgas being the 
gas pressure. Here, c r , c v and c z are Alfven sound speeds in three direction of cylindrical coordinate 
and define as 

where p mag z are the magnetic pressure in three directions. The angular transfer equation with 
considering outflow/wind and convection can be written as 

EiV-^-(mp) = --^-{Jvis) - --j-(Jcon) + rVT,c r -^(VY,c v ) + Ecvc^ - (^) (J) 
ar r ar r ar ar Zir ar 



where J vis and J con are viscous and convective angular momentum fluxes respectively which are 
defined as 

Jvis = -vYiT — (8) 
ar 

and 

Jcon = -V con ^+^U^-^) (9) 

ar 

Here, v is the kinematic viscosity coefficient, v = ac s H, with a being the constant Shakura & Sunyaev 
parameter, v con is the convective diffusion coefficient, and g is the index to determine the condition of 
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convective angular momentum transport. When g = 1, the flux of angular momentum due to convection 
is 

3 dQ 

J con = ~v con Tjr — (10) 
dr 

and when g = — 1/3, the convective angular momentum flux will be (Narayan et al. 2000) 

d(r 2 n) 

J con — V co11 YjT - (11) 

dr 

The last term on the right side of angular transfer equation, eq (|7]), represents angular momentum 
carried by the outflowing material(see e.g., Knigge 1999). 
By integrating over z of the hydrostatic balance, we have 

Q 2 K H 2 ~ -L Cr ^(VSc z )H - [c 2 + + cj)] = (12) 

By considering outward energy due to convection and energy loss by outflows, the energy equation 
will be ^ 

T,v r T— + - — {rF con ) = f(v + gv con )T,r 2 {—) 2 - -^rim w {r)v 2 K {r) (13) 
dr r dr dr I 

in above equation T is temperature, S is the specific entropy and F con is the convective energy flux 
which is defined as 

Fcon = -v con Y,T— (14) 
dr 

where 

T ds = j_d^^^dp (l5) 

dr 7 — 1 dr p dr 

here 7 is the radio of specific heats. As it mentioned, the last term on right hand side of energy equation 
is the energy loss due to wind or outflow (Knigge 1999). Depending on the energy loss mechanism, 
dimensionless parameter 77 may change. In our case, we consider it as a free parameter (Knigge 1999). 
Also in energy equation we still neglect the Joule heating rate. 

We adopt the assumptions of Narayan et al. (2000) and Lu et al. (2004) for the convective diffusion 
coefficient, v con , which is defined as 



M {-Nt ff y* (i6) 



L 2 

l^con ^ 

Here iV e // is the effective frequency of convective blobs and Lu is the characteristic mixing length. 
The effective frequency of convective blobs will be 

N 2 ff = N 2 + k 2 (17) 

with and k being the Brunt- Vaisala frequency and epicyclic frequency respectively, which are defined 

as 

p dr dr p 
amr 

Also, the characteristic mixing length Lu in terms of the pressure scale height, H p , can be written 

as 

L M = 2-V*l M H p (20) 

and 

H p = --¥— (21) 
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B v = —{v v B r - v r B v ), (25) 



where Im is the dimensionless mixing length parameter and its amount is estimated to be equal to V2 
in ADAFs (Narayan et al. 2000; Lu et al. 2004). The convective diffusion coefficient can also write in 
the form similar to usual viscosity of Shakura & Sunyaev (1973), 

l^con CX-conCs-B (22) 

where a con is a dimensionless coefficient that describes the strength of convective diffusion. The a con 
coefficient can be obtained by equations ([Tol l and (l22l 

<W=4ff7(-^//) 1/2 (23) 

Finally we can write the three components of induction equation, (B r , B v , B z ), to measure the magnetic 
field escaping rate, 

B r = 0, (24) 

d 

dr ' 

d , . v r B 7 

B z = - — (v r B z ) — . (26) 

dr r 

where B rtV)Z is the field scaping/creating rate due to magnetic instability or dynamo effect. Now we 
have a set of MHD equations that describe the structure of magnetized CDAFs. The solutions to these 
equations are strongly correlated to viscosity, convection, magnetic field strength (3 r>(PtZ and the degree 
of advection /. We seek a self-similar solution for the above equations. In the next section we will 
present self-similar solutions to these equations. 

3 SELF-SIMILAR SOLUTIONS 
3.1 Analysis 

We seek self-similar solutions of the above equations. Therefore, according to Narayan and Yi 1994, we 
can write similarity solutions as 

E(r) =c £ (-)^ (27) 
ro 



Vr (r) = -ciJ^(-)- L * (28) 
ro r 



1V(0=W — (-)"* (29) 
<?s(r)=ci(^)Q- 1 (30) 



c 



Pi CAT 

?, v ,,(r) = -f^ = 2/3 r ,,, z cl(™)(^)- 1 (31) 



r 



H{r) = c A r {—) (32) 
ro 

where constants cq, c\, C2, C3 and C4 will be determined later from the main MHD equation. £0 and 
ro are exploited to write the equations in non-dimensional forms and the constants P r ,tp,z measure the 
radio of the magnetic pressure in three direction to the gas pressure, i.e., /3 r)V)Z ~ Pmag r z /p ga s ■ 
In addition, the field scaping/creating rate B r ^^ z is assume to be a form of 

T .X / 11 \ 

B r ,ip,z = -E> r o iV o,zo( — ) 2 2 (33) 
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where £? r o,¥>o,zO are constant. 

By substituting the above self-similar solutions in the continuity, momentum, angular momentum, 
hydrostatic balance and energy equation of the disc, we obtain the following system of dimensionless 
equations to be solved for Cq, C\,c% C3 and c 4 : 



(34) 



\<?x = 4 - 1 - [(* - f ) + (s - I )A + (» + i)^]c; 



1 rh 
2>7 



C2 



3 1 1 

■ + 2^ tt + f«c n)C2C3C 4 + (s + 2) c 3(^^v) 1/2 



C4 = 



(s-^) 2 /3 r /3 z c2 + 4(l + /3 r +/3 v ) 



1/2 



/I 3\ f/ - 3 2I 9 2 / . 1 m 

(— ■ - + s - -) < (s - l)a C o„c 3 c 4 + cic 3 > = -/c 2 c 3 c 4 (a + gcw) - — S77 
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where m is the non-dimensional mass accretion rate which is defined as 

Mo 



If we solve the self-similar structure of the magnetic field escaping rate , we will have: 

= 



(35) 
(36) 

(37) 
(38) 

(39) 
(40) 



R 



1 



Oip 



7,GM, 



5/2 



c 3 



47rc S]o 
c 4 



■ 1 3 GM» 
So, = -(«- 2)2103^- 



47r^coSo 
c 4 



(41) 



(42) 



We can solve these simple equations numerically and clearly just physical solution can be inter- 
preted. Without mass outflow and magnetic field, i.e. s = I = rj = fi r = fi v = fi z = 0, the equations 
and their similarity solutions are reduced to the Narayan et. al solution (Narayan & Igumenshchiv & 
Abramowicz 2000). Also in the absence of outflow they are reduced to Zhang & Dai 2008. 

Now we can analysis behavior of the solutions in the presence of wind, convection and global mag- 
netic field. The parameters of our model are the standard viscose parameter a, the advection parameter 
/, the radio of the specific heats 7, the mass-loss parameter s, the degree of magnetic pressure to the gas 
pressure in three dimensions of cylindrical coordinate, (3 r , (3^ and (3 Z and /, 77 parameters corresponding 
to the rotation and non-rotation wind and energy loss by wind. 



3.2 Numerical Results 

In FigureQ] the surface density (co), the radial infall velocity (c 4 ), the rotational velocity (C2) and the 
sound speed (c 3 ) are shown as a function of toroidal magnetic field parameter /3 V for several values of 
wind parameter s, i.e., s = (dotted line, no wind), s = 0.1 (dashed line) and s = 0.2 (solid line). 
Four panels of figure[T]are set as a = 0.5, 7 = 1,, f3 r = (3 Z = 0.4, r) = I = 1 and / = 1. by adding 
the toroidal magnetic field parameter (3^, we see the surface density and rotational velocity of the disc 
increase, although the radial and sound speed both decrease. 
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On the other hand the radial flow decreases when the toroidal magnetic field becomes large. We can 
see that by adding the magnetic field influences (adding the rotation velocity will increase. This 
is because the disc should rotate faster than the case without the magnetic field which results in the 
magnetic tension. 

In figure[T]we also studied the effect of parameter s = on physical coefficients. As it mentioned, 
the value of s measures the strength of wind/outflow and a longer s denotes a stronger wind. We can 
see that for stronger outflow (s = 0), the reduction of surface density is more evidence. We see that the 
convective model of accretion flows with the presence of wind rotate more slowly than those without 
winds and wind leads to enhance accretion velocity. Also the sound speed of the disc increase for 
stronger outflows. 

Our variables as a function of magnetic field parameter (3 Z and several values of mass loses are 
shown in Figure [2] As it is seen, a strong z-component of magnetic field leads to increase both surface 
density and rotational velocity, while the radial infall velocity of materials and sound speed of the disc 
decrease. 

Figure [3] shows how the coefficients as depend on the magnetic field parameter in radial direction 
f3 r for several values of outflow parameter s. we see that the surface density and rotational velocity 
raise up when the magnetic field parameter f3 r increase. While the sound speed and accretion velocity 
will decrease. According to figures [TJ [2]and[3] by adding all components of magnetic field, the surface 
density and rotational velocity increase although the sound speed and radial infall velocity of the disc 
will decrease. 

In figure |4] we have plotted the convective parameter a con versus magnetic field parameters f3 r , f3 v 
and (3 Z for several values of s. As it is seen, when the magnetic field parameter j3 r become stronger, the 
convective parameter decrease (left panel). Also the toroidal component of magnetic field /L, decrease 
the a CO n too (middle panel), although by adding z-component of magnetic field the convective parameter 
a con increase (right panel). 

Figure [5] represent that the disc thickness enhance with toroidal magnetic field parameter (3 V and 
wind parameter s (middle panel). While the radial and vertical components of magnetic field decrease 
the vertical thickness of the accretion disc (left and right panels). 

4 CONCLUSIONS 

The CDAF model consistently represent radiatively inefficient accretion flows into black holes in the 
framework of viscose flow. We have considered the radial structure of RIAFs based on a self-similar 
analysis (CDAF). Some approximations were made in order to simplify the main equations. We assumed 
an axially symmetric, static disc with the a-prescription of viscosity, v = ac s H. A set of similarity 
solution was presented for such a configuration. We have extended Akizuki & Fukue 2006; Zhang 
& Dai 2008 and Faghei 2012 self similar solutions to present dynamical structure of the convection 
dominated accretion flows. We ignored the relativistic effects and the self-gravity of the discs. 

We have found out that by increasing all components of magnetic field , the surface density and 
rotational velocity increase although the sound speed and radial infall velocity of the disc decrease. Also 
we have shown that the existence of wind will lead to reduction of surface density as well as rotational 
velocity. Moreover the radial velocity, sound speed, advection parameter and the vertical thickness of 
the disc will increase when outflow becomes important in the radiation inefficiently accretion flow. 

In this manuscript we have follow the effect of Large-scale B-field on the structure convective this 
with wind and outflow, in the future it would be interesting to study how these effect would change the 
observational appearance of the flow. 

Although we have made some simplifying assumptions in order to treat the problem analytically, 
our self-similar solution shows Large-scale magnetic fields can really change typical behavior of the 
physical quantities of a hot accretion flow. Not only the surface density of the disc changes, but also 
the rotational and the radial velocities significantly change because of the magnetic fields. It means any 
realistic model for hot disc should consider the possible effects of the magnetic fields. 
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Fig. 1 Numerical coefficient as as a function of magnetic parameter j3 v for several values 
of s. The dotted, dashed and solid lines correspond to s = 0.0,0.1 and 0.2 respectively. 
Parameters are set as a = 0.5, 7 = 1, g = —1/3, f3 r = j3 z = 0.4, r) = I = 1 and / = 1, 
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rj = I = 1 and / = 1. 
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